Quantum spin dynamics of mode-squeezed Luttinger liquids in two-component atomic gases 
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We report on the observation of many-body spin dynamics of interacting, one-dimensional (ID) ultracold 
bosonic gases with two spin states. By controlling the non-linear atomic interactions close to a Feshbach reso- 
nance we are able to induce a phase diffusive many-body spin dynamics of the relative phase between the two 
components. We monitor this dynamical evolution by Ramsey interferometry, supplemented by a novel, many- 
body echo technique which unveils the role of quantum fluctuations in ID. We find that the time evolution of the 
system is well described by a Luttinger liquid initially prepared in a multimode squeezed state. Our approach 
allows us to probe the non-equilibrium evolution of one-dimensional many-body quantum systems. 
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Among the applications of ultracold atomic gases, atom 
interferometry stands out due to its potential for high preci- 
sion measurements Q. In atom interferometry, the physi- 
cal quantity of interest is measured in terms of the relative 
phase accumulated by the atomic wavefunction, subsequently 
mapped onto atomic populations for efficient read-out. Due 
to their intrinsic phase coherence and the possibility to cre- 
ate non-classical spin states for precision metrology, Bose- 
Einstein condensates (BEC) seem ideal candidates for such 
experiments. However, interatomic interactions mitigate this 
conclusion. For a two-component interacting BEC, it has been 
shown |2l using a single-mode approximation (SMA) that 
the relative phase between the two components undergoes a 
complicated evolution (Fig.[T]c-e), creating quantum correla- 
tions |4| while single -particle coherence is suppressed. There- 
fore, this dynamics is often termed phase diffusion. 

In this work, we investigate such an interaction induced dy- 
namics in quasi- ID two-component quantum gases by moni- 
toring the loss of coherence in a Ramsey-type interferometer 
sequence. In order to distinguish different contributions af- 
fecting the coherence through the spin or spatial wave func- 
tions, we employ a novel many-body spin echo sequence us- 
ing a Feshbach resonance to adjust sign and magnitude of the 
atomic interactions. When applied to a single spatial mode 
BEC, this spin echo would lead to full revivals of coherence, 
which are not observed in our experiment. In contrast, quan- 
tum fluctuations play a key role for ID interacting systems 
J5] [6), which must necessarily be described as multimode 
quantum gases, as during the dynamical evolution higher en- 
ergy modes become populated [7 1. The Luttinger liquid (LL) 
formalism |8 9 1, which reduces the interacting problem to an 
effective low-energy model of decoupled harmonic oscillator 
modes, provides such a description. We show theoretically 
that our preparation sequence amounts to producing a multi- 
mode squeezed state in the spin excitation modes of the LL os- 
cillators, with each oscillator itself prepared in a well-defined 
mode-squeezed state, and remaining in a squeezed state at all 
times [ 1 1 . Monitoring the phase dynamics of this strongly 
non-equilibrium state allows to probe fundamental aspects of 
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FIG. 1: (a) Array of quasi- ID spinor systems. (b,c) Within each 
tube, a coherent spin-state is created exhibiting Gaussian distributed 
fluctuations of the mean spin, (c-e) Time evolution of the initial CSS 
under non-linear interactions re-distributing the initial Gaussian fluc- 
tuations towards increased phase fluctuations assuming a SMA. 

ID physics, namely the competing dynamics of the (quasi-) 
condensate fraction (zero momentum mode) and of the low- 
energy excitations, highly relevant for squeezing experiments 
in ID configurations [IT). From our model we find that only 
the lowest oscillator mode shows the familiar revival dynam- 
ics, whereas the full model leads to the partial revivals that we 
observe experimentally. 

The system we consider is an array of two-component 87 Rb 
spinor gases confined to quasi-lD traps (tubes). We experi- 
mentally realize this system by loading a 87 Rb BEC of around 
2.8 x 10 5 atoms into a 2D-optical lattice [12|, hence creating 
a 2D array of ID degenerate quantum gases (see Fig.[T^). The 
lattice laser wavelength is A — 843 nm, and the radial and ax- 
ial trap frequencies within each tubes are a> r ~ 2n x 42 kHz, 
and bi ax as 2n x 90 Hz, respectively, from which we calculate 
a mean number per tube of N - 60. 

In order to extract information about the phase dynam- 
ics, we experimentally monitor the coherence of the system 
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by recording interference fringes in a Ramsey-type interfer- 
ometer. In the following we use the well known analogy 
with a quasi spin- 1/2 system in order to describe our two- 
component gas. Starting from a spin-polarized ensemble in 
state ||) = \F = \,m.f = +1), we use a two-photon 7r/2-pulse 
combining a microwave and a radio frequency photon to cou- 
ple this state to the | f) = |2, -1) spin state and bring each atom 
into the single-particle superposition (| T) + 1 I))/ V2. This pre- 
pares a coherent spin state (CSS) within each tube with ex- 
pectation value of the magnetization (m z ) = and variance 
(m 2 ) - N/2, with m z - - «| (see Fig. [IJ), cf. T\3\ ). In 
order to observe interaction driven effects, we let the system 
evolve for a given time at a particular value of the inter-spin- 
state interaction strength, selected by using a Feshbach reso- 
nance around B — 9.12 G. Thereby the inter-species scattering 
length an can be changed by a few 10% from its background 
value lTT4l [151 . where ay is the s- wave scattering length for 
collisions between atoms in spin states i and j. 

After this time evolution, a final 7r/2-pulse with phase 6 rel- 
ative to the first pulse is applied, mapping the final relative 
phase onto populations of spin states | f) and | J,) that are read 
out using state-selective absorption imaging. In the absence 
of interactions and dephasing, such a sequence results in si- 
nusoidal Ramsey fringes in the relative population Nf /N to t as 
a function of 6. Experimentally the coherence is quantified 
through the visibility of the Ramsey fringe 

^L = I(l+<y(,)xco S (0)), (1) 

which is used to fit the experimental data and extract 'V(f) for 
a specific interaction time. Far from the Feshbach resonance 
(B = 8.7 G), where the effect of the phase dispersion can be 
neglected, we measure a e -1 -decay time fd ec = 54 ms, which 
can be attributed to residual single particle decoherence ef- 
fects, e.g. caused by magnetic field fluctuations. 

Close to the Feshbach resonance, however, we find a 
markedly faster decay of the Ramsey contrast. In Fig. [2] we 
monitor such a behavior of the Ramsey fringe over time for 
two magnetic fields located almost symmetrically around the 
center of the Feshbach resonance. Such a behavior can be 
expected from enhanced phase diffusion due to increased in- 
teractions near the resonance. Phase diffusion results from a 
spread in the distribution of populations which are converted 
into phase fluctuations by the non-linear interactions during 
the evolution, see Fig. [TJ;-e. In the simplest case where all 
atoms occupy the same orbital wave function, the Ramsey 
fringe contrast decays according to 

^sMA(f)^expf-^ 2 <m2) f 2 J. (2) 

For the initial state we prepare, the population variance 
(m 2 ) - N /2 leads to a phase uncertainty (A0o) 2 = 2/N * 
0.033 of the collective spin vector in the equatorial plane 
(see Fig. 1) in each tube, and a phase spreading time scale 



h ~ 1 / Of yN)- The parameter^, related to the second deriva- 
tive of the chemical potential [3|, is directly proportional to 
the difference a s - (a^ + - 2an)/2. Far from the reso- 
nance, all three scattering lengths a^, a^, and a-|-j_ are approx- 
imately equal, so that y ~ and interaction-induced phase 
spreading can be neglected. However, near the Feshbach 
resonance, the change of inter-species scattering length can 
lead to a significant non-linear interaction energy. Following 
Refs. HUEIQI], we estimate X ~ 2ttx4.6Hz for B = 9.131 G 
and our trapping parameters with an atom number of N ~ 60. 
Although this value, together with the observed decoherence 
rate, is roughly on the order of the observed rate at which the 
coherence is lost in the quasi- ID regime investigated here, the 
pure SMA Eq. Q cannot explain our experimental observa- 
tion in Fig. [2] Close to the resonance we lose up to 50% of 
the atoms due to inelastic collisions. However, as these col- 
lisions usually remove atoms from both spin states symmet- 
rically, they do not modify the magnetization of the system 
and thus only weakly influence the dynamical evolution of the 
coherence for our measurement times ||3] Q~8] . 

In contrast to the simple model of Eq. Q which predicts 
a symmetric decay around the resonance, we systematically 
observe a faster drop of contrast below the resonance. In fact, 
changing the sign of the effective interaction strength has se- 
vere consequences on the spatial wave function of the atoms. 
Below resonance, the inter-species repulsion is stronger than 
the intra-species repulsion (y < 0), and the system becomes 
dynamically unstable towards demixing of the two species. 
This reduces the Ramsey fringes visibility below the reso- 
nance, which cannot be distinguished from the effect of the 
coherent phase diffusion dynamics. 

In order to separate the effects of phase diffusion from other 
mechanisms reducing the Ramsey fringe contrast, we apply a 
many-body spin echo operation after an initial evolution time 
T, similar to the one used in cavity quantum electro dynamics 
experiments 1 19 1. We stress that our echo technique is acting 
on the many-body quantum state, thereby extending previous 
theoretical work on echo operations neglecting phase diffu- 
sion ll20l . Such a many -body spin echo operation is performed 
by first holding the sample for a time T = 6 ms at a magnetic 
field B\ = 9.131 G above the Feshbach resonance (y > 0), 
and subsequently jumping below the resonance (x < 0). This 
operation effectively changes the sign of the non-linear inter- 
action parameter^, while heating or atom loss can be avoided 
lfT3l . In a SMA one would expect this sequence to correspond 
to a perfect time reversal, leading to a full revival of the con- 
trast after another interaction time T according to 

TWO = ^ exp(-^<m 2 )(?-2r) 2 J. (3) 

This contradicts our observation of only partial revivals, 
shown for two spin echo sequences in Fig. [3] 

In order to explain our observation, we model our system in 
a LL approach going beyond the usual SMA J2l[7|. A drastic 
simplification follows from the near equality a\\ ~ flff, which 
results in a decoupling of elementary excitations into almost 
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FIG. 2: Ramsey fringe contrast drop for a time evolution at at 
B = 9.106G (O) and B = 9.131 G (•). The dashed line indicates 
the independently measured decoherence far away from resonance 
(fdcc = 54 ms). The solid lines are predictions of our LL model with 
the phase width (A(p LL ) 2 = (A0 O ) 2 ~ 0.033 fixed. The dashed dotted 
line is a prediction of Eq. |2]( based on SMA. For values see text. 



independent density and spin fluctuations. The latter can be 
described in terms of two conjugate fields, m z and $> s , describ- 
ing respectively fluctuations of the local magnetization and of 
the relative phase. At long-wavelengths, the spin-part of the 
Hamiltonian reads 



dx 



g s m z 



4M 



(4) 



with n tot = rif + n± the linear density. For a uniform ID 
system of length L, the fields m, and cf> s can be expanded 
in terms of momentum eigenmodes a q ,at as lETTl (f> s (x) = 
(po + T jq (2qLKIn)- ll2 e-^ 12 ^ sgn(^) \e iqx a q + h.cl Each mode 
is characterized by the wave vector q and frequency u> q (t) - 
v s (t)\q\, where v s (t) = (g s (t)n tot /M) 1 ^ 2 is the spin velocity and 
g s (t) denotes the spin coupling constant. The sum over LL 
modes exclude the zero mode, and are restricted to values of 
q below a cut-off momentum q c ~ £ h where is the healing 
length. The interactions are encoded in the LL parameter K 
which can take values ranging from 1 (the so-called Tonks- 
Girardeau limit lfT71l22lD to infinity (noninteracting gas). For 
weakly interacting bosons [9] K * n/[y/y(l - yfylln) 1 ^ 2 ] 
where y - 2a s Mcj r /hn tot . In our experiments lfT31 y ~ 
0.1 - 0.2 and K ~ 5 - 8 for different data sets. The con- 
tribution of the zero energy mode <f>o is identical to that in 
the SMA, as described before [2]. The dynamics of the low- 
energy LL excitations on top of the zero mode dynamics is 
that of a collection of independent harmonic oscillators with 
(time-dependent) frequencies cj q . 

In order to use the LL decomposition to compute the time 
evolution, we need to identify how to describe the initial state 
in terms of those LL modes. Our experimental scheme ideally 
corresponds to an instantaneous projection of the spin state 
(initially polarized in | J,)) onto a state with zero relative phase 
directly after the first 7r/2-pulse, (j> s (x)\if/(0)) w 0. The connec- 
tion with the LL formalism is done by identifying the initial 
CSS as a multimode squeezed state for the elementary spin 



excitations, 

|-A(0)> = £(KD|0> = ]~[ ^/(l -| W? | 2 )exp( W? 4at ? )|0> (5) 

i 

where S is a squeezing operator, with the factor w q = (1 - 
<Xq)/(l + o: q ). Here, a q = A(f>u_\q\/\q c \ is a mode squeezing 
parameter, and the phase variance (Atyn) 2 = (A0o) 2 = %■ 
In reality, due to various experimental imperfections (e.g. un- 
known temperature), the exact width of the prepared squeezed 
state is unknown and cannot be determined independently. We 
still consider the initial state as a squeezed state of the LL os- 
cillators, with a fitted Atpn of the squeezed state. 

The time evolution of a squeezed state under the LL Hamil- 



tonian amounts to the replacement a q 



j exp(-icj q t). In 



addition, the reversal of the sign of interaction at time t — T 
amounts to a sign reversal of the spring constant of each 
LL harmonic oscillator. We are able to compute this time 
evolution exactly iTPJI . using the formalism of harmonic os- 
cillators with time-dependent frequencies co q (t) f23j. Here, 
we concentrate on the comparison between the predictions 
of the calculations and the experimental results. From our 
model we are able to calculate the coherence factor < V(f) = 

^ e \l Jo dx(t//(t)\e'^'\i//(t))^j, measuring the relative phase S 
between ||) and |J,). In the LL formalism, the coherence fac- 
tor can be generally written as 'V(f) = *Vsma(0 X *Vq*o(t), 



where 'Vsma(0 is given by Eqs. (2j3 1, and the term describ- 
ing the contribution of the q + modes to the decay of con- 
trast < Vq?o(t) is known explicitly Ifl3l . The combination of 
these effects leads to the typical behavior of 'V(f) illustrated 
in Fig. [3] where we quantitatively compare the experimental 
results with the predictions of our LL model. The interaction 
parameters of the LL model were determined from the micro- 
scopic data 03], and in our computations we only consider 
the density in the central tubes computed in the Thomas -Fermi 
approximation lfl8l . Overall, we find very good agreement be- 
tween the LL model and the experimental data using A(pn as 
the only fit parameter. For our data we find values that are of 
the same order of magnitude as the initial width. For longer 
times (f > 20 ms), the model deviates from the measured 
data. This breakdown is due to the phenomenon of demix- 
ing discussed above, when the excitations become so strong 
that density fluctuations are significant [6|. Its timescale can 
be estimated as the time required for the formation of random 
magnetization domains, when J^d^l 2 ) ~ N. For our experi- 
mental parameters, it is of the order 20 - 25 ms which is of 
the order of 1 /%. 

From this analysis we find indeed that although the zero 
mode evolution is perfectly refocused by changing the sign of 
g s , the non-zero modes still undergo dephasing even under the 
echo sequence. The reason for this is the kinetic energy term 
in Eq. Q, unaffected by the spin echo. Hence, the reversal is 
exact for the q — mode, significant for low-lying spin waves 
with q ~ l/L, but increasingly less efficient for higher lying 
spin waves modes with L 1 <K q < q c . Note that a full revival 
could be in principle achieved by also reverting the second, 
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FIG. 3: Ramsey fringe visibility versus time for holding above the 
Feshbach resonance, B = 9.131 G (•), and with time reversal (spin 
echo) after 6.7 ms (O) jumping from B = 9.131 G to (a) B = 9.101 G 
and (b) B = 9.090 G. The solid and dashed lines are calculation from 
our model with (A0 Li ) 2 = 0.04 (a), and (A0 Li ) 2 = 0.025 (b). For 
B = 9.131 G, a s = 0.17 a-pj and we compute / = 0.165 and A* = 8.0. 



kinetic energy term of the Hamiltonian in Eq. Q. This could 
be realized by inducing a negative effective mass, e.g. through 
a weak optical lattice along the direction of the tubes Il24"ll . 

In conclusion, we have studied the phase dynamics of 
quasi-lD two-component quantum gases with adjustable in- 
teraction. For strong interactions we observe an accelerated 
decay of coherence in the system. By application of a novel 
many-body spin echo technique we are able to reverse the in- 
teraction driven dynamics, leading to a partial revival of the 
coherence in the system. We attribute this revival to the dy- 
namics of the ground state mode, reflecting the coherent na- 
ture of the phase diffusion dynamics. This is supported by 
quantitative comparison to our LL model. The missing frac- 
tion of coherence in the revival is also quantitatively in agree- 
ment with our model and demonstrates the importance of the 
dynamical evolution of higher lying modes in ID systems. 
Our experiment shows that quantum fluctuations are a cru- 
cial component in the discussion of phase diffusion dynamics 
[2 1 and spin-squeezing Bl fTTTl in low dimensional systems. 
While our work demonstrates that these quantum fluctuations 
fundamentally limit the performance of atom interferometers 
in ID, it also indicates an avenue to overcome such limitations 
by inverting both interaction and kinetic energy terms simul- 
taneously during the interferometer sequence. 
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